1/(N — 1) expansion based on a perturbation theory in U 
for the Anderson model with TV-fold degeneracy 
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We study low-energy properties of the TV-fold degenerate Anderson model. Using a scaling that 
takes u = (TV — 1)U as an independent variable in place of the Coulomb interaction U, the perturba- 
tion series in U is reorganized as an expansion in powers of 1/(TV — 1). We calculate the renormalizcd 
parameters, which characterize the Kondo state, to the next leading order in the 1/(TV — 1) expansion 
at half-filling. The results, especially the Wilson ratio, agree very closely with the exact numerical 
renormalization group results at TV = 4. This ensures the applicability of our approach to TV > 4, 
and we present highly reliable results for nonequilibrium Kondo transport through a quantum dot. 
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The Anderson impurity has been studied extensively as 
a model for strongly correlated electrons in dilute mag- 
netic alloys, quantum dots, and also for bulk systems 
in conjunction with the dynamical mean- field theory 
For quantum dots, the nonequilibrium Kondo effect can 
occur when a bias voltage is applied between two leads. 
A universal Fermi-liquid behavior [H-Q has been closely 
examined at low energies for the steady current 0- 11 
and shot noise fl2- 17|. 

Orbital degeneracy in the impurity states also affects 
the nonequilibrium properties at low energies. Recently, 
Mora et al. lfj have succeeded to express the current 
noise in terms of the Fermi-liquid parameters [2- 5] in an 
SU(iV) Kondo regime, where the Coulomb repulsion U 
is so large that charge fluctuations arc suppressed near 
the impurity with TV-fold degeneracy. A complemen- 
tal expression that takes into account the fluctuations 
at half-filling has been presented in our previous work 
18[. In this case, the corrections due to finite U en- 
ter through the Wilson ratio R, which is a correlation 
function defined with respect to the equilibrium ground 
state, and through the width of the Kondo resonance 
A. Therefore, explicit values of these two parameters, 
R and A, are required to study the low-energy trans- 
port thoroughly. The exact numerical renormalization 
group (NRG) approach is still applicable to multi-orbital 
systems. It practic ally works, however, for small degen- 
eracies TV < 4 [3, [l9[ , which for TV = 2 corresponds to 
the spin degeneracy. Therefore, alternative approaches 
are needed to explore the large degeneracies at TV > 4. 

In this Letter, we propose a systematic approach to 
calculate correlation functions at TV > 4, using a scaling 
that takes u = (TV — 1)U as an independent variable in 
place of U. Here, the factor TV — 1 corresponds to the 
number of different impurity states, with which a local 
electron in the impurity site can interact. With this scal- 
ing, the perturbation series in U can be reorganized as an 
expansion in powers of 1/(TV — 1), using a diagrammatic 



classification similar to the one for the TV-component tp 4 
model (ioj . However, our approach is completely dif- 
ferent from the usual 1/TV expansion and non-crossing 
approximation, which are constructed on the basis of the 
perturbation expansion in the hybridization matrix ele- 
ment v v [21 - 23 1 . We calculate R and A up to the next 
leading order terms in the 1/(TV — 1) expansion at half- 
filling, and find that the results agree very closely with 
the NRG results at TV = 4, where TV is still not so large. 
Particularly, the Wilson ratio shows an excellent agree- 
ment over the whole range of U. The early convergence 
of the expansion implies that our scaling procedure effi- 
ciently captures the orbital effects, and ensures the ap- 
plicability to TV > 4. This enables us to present highly 
reliable results for the nonequilibrium steady current and 
shot noise for TV > 4. Our approach could have wide ap- 
plication to quantum impurities, and could be used as a 
solver for the dynamical mean-field theory 24 1. 

The Hamiltonian for the TV-fold degenerate Anderson 
model connected to two leads (y = L, R) is given by 

H = Ho+Hu, Hu = - Un d m n dm>, (!) 

— ^ ] ^ ] / ^ e 6 c \um C cwm + ^ ] e d C ^ln < ^m^ 
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v=L,R m=l 

Here, S m creates an electron with energy €d in orbital m 
at the impurity site, n dm — <v m d m , and m (= 1, 2, ■ • • , TV) 
includes the spin degrees of freedom. c\ vm creates a con- 
duction electron with energy e and orbital m in lead v, 



and is normalized as {c e 

The linear combination ip um = J D D de^/pc euml with 
p = 1/(21?), couples to the impurity level via the hy- 
bridization matrix element v , and A = Tl + with 



2 



= Trpv 2 . We consider the parameter region where A, 
Cd, and U are much smaller than the half band width D. 

We use the imaginary- frequency Green's function that 
takes the form G(iuj) = [iu — e d + iA sgnw — Y,(iuj)}~ 1 
for | a; | <C D. The behavior of the self-energy E(iw) 
for small w determines the enhancement factor for the 
linear specific heat 7 = 1 — dT,(icj) /d(iui)\ u =o, and the 
renormalized parameters z = 1/7, ed = z[ed + S(0)], 
and A = zA. The average number of local electrons 
can be deduced from the phase shift 6 = cot -1 (e d /A), 
using the Friedel sum rule, (ridm) — 5 /it. The enhance- 
ment factor for the spin susceptibility and that for the 



charge can be written in the form x s 



X 77 



Xmm' and 



Xc = Xmm + (N-1) Xmm' for m ^ ™.' . These susceptibil- 
ities can be deduced from the self-energy and four-point 
vertex function T mm ,. m , m (iuj 1 , iw 2 ; «w 3 , viw 4 ) for m =^ m' , 
using the Ward identities [5(, 



sin 2 S 
ttA 



- mm ;m m 



, m (0,0;0,0). (3) 



Furthermore, U = z 2 T 



mm ;m m 



(0,0; 0,0) corresponds to 



the residual interaction between the quasi-particles. 

The Wilson ratio R parameterizes how far the system 
is away from the Kondo limit, and plays a central role 
for finite U, 



R=^ = l 
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^ = l-gsin 2 S. (4) 



Here, the scaling factor N — 1 is introduced to the renor- 
malized interaction U and the bare one U, such that 



9=(N-1) 



U 
ttA 



(N-l) 



U 
ttA 



(5) 



In the following we consider the particle-hole symmetric 
case, where e d = —(N — l)U/2 and 5 = tt/2. In this 
case, the renormalized coupling takes a value in the range 
<g < 1. It approaches to j? — > 1 in the limit of g — > 00 
as the charge fluctuation is suppressed \ c — > 0. 

We calculate 7 and r mm/ . m/m (0, 0; 0, 0) perturbativcly 
to order U 3 and U 4 , respectively, by extending Yamada's 
calculations for TV = 2 0] to general N [l8j|, and obtain 
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(6) 
(7) 



Here, C( x ) i s the Riemann zeta function, which disap- 
pears at N = 2 where the impurity has only the spin de- 
generacy For TV > 2, 7 and g are no longer even nor 
odd function of U. We see in Eqs. (j6]) and that the 
coefficients in the perturbation series can be expanded 
in powers of 1/(N — 1). Thus, the perturbation series 
in g can be reorganized as an expansion with respect to 
1/(N — 1). If the N — > 00 limit is taken at fixed g, 
then the right hand side of Eq. (JB|) approaches to an al- 
ternating geometric series in g, and 7 approaches to the 
noninteracting value 7 —J- 1. We will see later that these 
are true for all order in g, and the asymptotic forms of 
Eqs. ([6]) and ((7]) in the large N limit are given by 



9 



1+9 



O 



N - 1 



7=1 + 



N — 1 



(8) 



The corrections due to finite N can be extracted, using 
a diagrammatic representation of the perturbation in U . 

The leading order contributions in the 1/(N— 1) expan- 
sion arise form a series of the bubble diagrams indicated 



in Fig. [1] and the sum of these diagrams corresponds to 

(9) 



4>(iu) gTrAn(icj) / 1 
^bubM = — — - + — — — + O 
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Here, Xo( iuJ ) = - J^rG (iuj + icj')G n (iui'), and 
G (iuj) = [iu - E d + liAsgno;] -1 with E d = [25|. Thus 
Xo( iuJ ) = ^ 2 \x^+\x\) with x = uj/A. The propagator 
U huh (iio) contains not only the leading order, but also 
higher order contributions in the 1/(N — 1) expansion. 
This is because the orbital indices for adjacent bubbles 
have to be different, and summations over internal m's 
arc not independent. The order 1/(N — 1) contributions 
to the vertex and self-energy come from the diagrams 
shown in Fig. [2] 

To calculate the renormalized coupling constant g to 
order 1/(N - 1), we need r mm ,. m , m (0, 0; 0, 0) to order 
1 / (N —l) 2 as g has a scaling factor N — 1 defined in Eq. 



FIG. 1. The leading order diagrams in the 1/(N — 1) expan- 
sion. The wavy and solid lines indicate the Coulomb repulsion 
U and unperturbed Green's function Go, respectively. The 
double wavy line represents the sum of the bubble diagrams, 
and corresponds to Wbub(ii^) given in Eq. ((9J. 





FIG. 2. The diagrams which provide the order 1/(JV — 1) 
contributions with some higher order corrections [see Eq. (JSJ] . 



((SJ). The order 1/(N — l) 2 contributions to the vertex 
function arise from the diagrams shown in Fig. [3J and 
from the order l/(N — l) 2 component of the vertex di- 
agram in Fig. [2J Summing up all these contributions, ~g 
can be expressed in the form that is exact up to terms of 
order 1/{N - 1), 
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Here,Z ( 5 ) = wAf^ {G M} 2 II(tw), and IV' = 7V-1. 
This formula shows the correct asymptotic form in both 
the weak and the strong coupling limits: g — g for g — > 0, 
and <jr — > 1 for g — >• oo. Thus, Eq. (fTTT) can also be 
regarded as an interpolation formula for the Wilson ratio 
as -R-l = g/(N-l) at half-filling. The order 1/(N-1) 
results for g show an excellent agreement with the NRG 
results for N — 4 as indicated in Fig. [5] (a). 

To obtain Eq. , the parameter 7 in the denomina- 
tor has been taken into account up to order 1/(N — 1), 
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N — 1 
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We also calculate, 7 



the order l/(iV — l) 2 contri- 



butions which arise from the diagrams shown in Fig. 0] 
and from the higher order component of the self-energy 
diagram in Fig. [2] 

Figure [5] (a) shows a comparison between the NRG 
O [ill and the 1/(JV - 1) expansion results for N = 
4. We see the very close agreement, especially for g. 
Although the order 1 / (N — 1) results are slightly smaller 
than the NRG results, the two curves for g almost overlap 
each other over the whole range of g. The deviation must 
decrease as N increases. Therefore, the order 1/(N — 1) 
formula for g given in Eq. (fTTj) provides almost exact 
numerical values for N > 4. We also see in Fig. [S] (b) the 
value that 5 can take is bounded in a very narrow region 




FIG. 3. The order 1/(N — 1) diagrams for the vertex function 



mm :m m 



(0,0; 0,0) for m^m'. 




FIG. 4. The order 1/(N - l) 2 self-energy diagrams which 
contribute to the renormalization factor z (= I/7). 



between the curve for N = 4 and that for the N — > 00 
limit. As TV increases, g varies rapidly towards the value 
for the large N limit. The order l/(N — l) 2 results for 
the renormalization factor z, shown in Fig. [S] (a), also 
agree with the NRG results for A = 4at<?<3.0, or 
cquivalently g < 0.8, from the weak to the intermediate 
coupling region where g is still not converged to 1.0, the 
value for the strong coupling limit. Therefore, away from 
the strong coupling regime the Kondo energy scale, A = 
zA, can be deduced reasonably from the order 1/(N — l) 2 
results. 

The 1/(N — 1) expansion can be applied fruitfully to 
nonequilibrium transport at finite U. To be specific, we 
choose the lead-dot couplings and chemical potentials to 
be symmetric: Tl = Tr and fj,£ = —j-iR (= eV/2). In this 
case, an exact expression can be derived for the retarded 
Green's function at low energies up to order w 2 , T 2 , and 
(eV) 2 gGl, 



uj + iA + i 



2(iV-l)A 



uj2 + 3( eV y + (irTy 



(13) 



The differential conductance for the current through the 
impurity can be deduced from G r (uj), using the formula 
by Meir- Wingreen (2(| and Hershfield 27 1 , 
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The low-energy behavior is characterized by the two pa- 
rameters, g in the coefficients and A the energy scale, 
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FIG. 5. (Color online) (a): <j and z versus g for TV = 4. The 
curve with the circles represents the NRG results. The red 
dotted line represents the order 1/(TV — 1) results for fl, and 
the blue dashed line the order l/(N — 1 
g vs g for TV = 2 (Bethe ansats Q), TV 
(order 1/(TV — 1)), and for TV — ► oo where p — > g/(l + g) 



2 results for z. (b): 
4 (NRG), TV = 6 



which depend on TV. Figure [5] (a) shows the ratio of 

c v to c T as a function of g for several TV, using Eq. 

(fTTI) for TV > 6. The ratio takes a value in the range 

3/4 < c v /c T < (3/4) (AT + 4)/(TV + 1) The order 

1/(TV — 1) results for lj are numerically almost exact for 

TV > 4 as mentioned, and thus the results shown in Fig. 

|5] capture orbital effects correctly. 

As another application of Eq. (jTTj) . we also consider 

the shot noise S = Jdt (6J(t)SJ(0) + SJ(0)SJ(t)), where 

5J(t) = J(t) — (J) is the current operator. At T = 0, 

S has been calculated to order (eV) 3 for the symmetric 

Anderson model for TV = 2 0, an d f° r general TV: 

Q — Ne2 1 
h 6 ' 



1 



eV m. The Fano factor F b 



9T \ (eV_ s 
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is defined as the ratio of S to the backscattering current 
J;, = NeV/h— J, and has been obtained in the form [IF 
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(16) 



It takes a value in the range 1 < F b < (N+8)/(N+4). In 
Fig. E](b), the order 1/(TV — 1) results for F b are plotted 
as functions of a fo r TV > 6, together with the exact re- 
sults for TV < 4 [18|]. As TV increases, g converges rapidly 
to the value, g ~ g/(l +g), for the large TV limit, as men- 
tioned in the above. Thus, for TV > 8, the TV dependence 
is determined essentially by the factor 1/(TV — 1), seen 
explicitly in Eq. ([TBI . The l/(Af — 1) expansion can also 
be applied to the full counting statistics [29| . 

In conclusion, we have described the 1/(TV— 1) expan- 
sion approach based on the scaling defined in Eq. (|5|). 
The next leading order results for g, which at half-filling 
corresponds to g = (TV — 1)(R — 1), can be expressed in 
the form of Eq. (|TT|) . We find that this formula inter- 
polates almost exactly between the weak and the strong 
coupling limits for TV > 4. The 1/(TV — 1) expansion 
can be extended to explore the particle-hole asymmet- 
251 ] - Furthermore, it provides a well-defined and 



ric case 



controlled way to take into account the fluctuations near 
the TV — > oo fixed point of many fcrmion systems with 
two-body interactions. 
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